*Stata .do file for "When is 'Time's Up'? The Influence of Severity and Costs/Benefits on Perceptions of Whistleblowing"

****
****


*front-end of paper text references

	*get culled data set
clear
use "/.../WB replication data.dta"

	*# of men and women in experiment
tab gender

	*N's and %'s of DK reponses
		*sexual harassment, female target
tab MIA320
		*harassment, male target (note: post-election)
tab MIA405
		*sexual harassmwent bystander
tab MIA321
		*chemicals in water
tab MIA322
		*embezzling, bystander, $500-5k
tab MIA323
		*embezzling, bystander, $1k-50k (note: post-election)
tab MIA407
		*embezzling, bystander, $2k-100k (note: post-election)
tab MIA408
		*leaking
tab MIA324
		*exam cheating, rewards complaintant
tab MIA325
		*exam cheating, risk of losing job if reported
tab MIA326
		*supervisor threatens public interest, job security protected by law
tab MIA327
		*supervisor violates federal law, positive reward for complaintant and issue will be resolved
tab MIA328

	*n's by treatment condition
		*sexual harassment, female target
tab MIA320_treat
		*sexual harassmwent, female target, bystander
tab MIA321_treat
		*harassment, male target (note: post-election)
tab MIA405_treat
		*chemicals in water
tab MIA322_treat
		*embezzling, bystander, $500-5k
tab MIA323_treat
		*embezzling, bystander, $1k-50k
tab MIA407_treat
		*embezzling, bystander, $2k-100k
tab MIA408_treat
		*leaking
tab MIA324_treat
		*exam cheating 1
tab MIA325_treat
		*exam cheating 2
tab MIA326_treat
		*supervisor threatens public interest 1
tab MIA327_treat
		*supervisor threatens public interest 2
tab MIA328_treat


****
****


*main findings: between subjects effects

	*get culled data set
clear
use "/Users/caseyklofstad/Dropbox/2018 CCES UM/paper 1 w-blowing and metoo/WB replication data.dta"

	*clean DVs: flip to run def not - def yes; put "don't know/not sure" in the middle
		*sexual harassment, female target
generate MIA320_flip=.
replace MIA320_flip=5 if MIA320==1
replace MIA320_flip=4 if MIA320==2
replace MIA320_flip=2 if MIA320==3
replace MIA320_flip=1 if MIA320==4
replace MIA320_flip=3 if MIA320==5
tab MIA320 MIA320_flip
		*sexual harassment, female target bystander
generate MIA321_flip=.
replace MIA321_flip=5 if MIA321==1
replace MIA321_flip=4 if MIA321==2
replace MIA321_flip=2 if MIA321==3
replace MIA321_flip=1 if MIA321==4
replace MIA321_flip=3 if MIA321==5
tab MIA321 MIA321_flip
		*sexual harassment, male target
generate MIA405_flip=.
replace MIA405_flip=5 if MIA405==1
replace MIA405_flip=4 if MIA405==2
replace MIA405_flip=2 if MIA405==3
replace MIA405_flip=1 if MIA405==4
replace MIA405_flip=3 if MIA405==5
tab MIA405 MIA405_flip
		*chemicals in water
generate MIA322_flip=.
replace MIA322_flip=5 if MIA322==1
replace MIA322_flip=4 if MIA322==2
replace MIA322_flip=2 if MIA322==3
replace MIA322_flip=1 if MIA322==4
replace MIA322_flip=3 if MIA322==5
tab MIA322 MIA322_flip
		*embezzling, bystander, $500-5k
generate MIA323_flip=.
replace MIA323_flip=5 if MIA323==1
replace MIA323_flip=4 if MIA323==2
replace MIA323_flip=2 if MIA323==3
replace MIA323_flip=1 if MIA323==4
replace MIA323_flip=3 if MIA323==5
tab MIA323 MIA323_flip
		*embezzling, bystander, $1k-50k
generate MIA407_flip=.
replace MIA407_flip=5 if MIA407==1
replace MIA407_flip=4 if MIA407==2
replace MIA407_flip=2 if MIA407==3
replace MIA407_flip=1 if MIA407==4
replace MIA407_flip=3 if MIA407==5
tab MIA407 MIA407_flip
		*embezzling, bystander, $2k-100k
generate MIA408_flip=.
replace MIA408_flip=5 if MIA408==1
replace MIA408_flip=4 if MIA408==2
replace MIA408_flip=2 if MIA408==3
replace MIA408_flip=1 if MIA408==4
replace MIA408_flip=3 if MIA408==5
tab MIA408 MIA408_flip
		*leaking
generate MIA324_flip=.
replace MIA324_flip=5 if MIA324==1
replace MIA324_flip=4 if MIA324==2
replace MIA324_flip=2 if MIA324==3
replace MIA324_flip=1 if MIA324==4
replace MIA324_flip=3 if MIA324==5
tab MIA324 MIA324_flip
		*exam cheating, rewards complaintant
generate MIA325_flip=.
replace MIA325_flip=5 if MIA325==1
replace MIA325_flip=4 if MIA325==2
replace MIA325_flip=2 if MIA325==3
replace MIA325_flip=1 if MIA325==4
replace MIA325_flip=3 if MIA325==5
tab MIA325 MIA325_flip
		*exam cheating, risk of losing job if reported
generate MIA326_flip=.
replace MIA326_flip=5 if MIA326==1
replace MIA326_flip=4 if MIA326==2
replace MIA326_flip=2 if MIA326==3
replace MIA326_flip=1 if MIA326==4
replace MIA326_flip=3 if MIA326==5
tab MIA326 MIA326_flip
		*supervisor threatens public interest, job security protected by law
generate MIA327_flip=.
replace MIA327_flip=5 if MIA327==1
replace MIA327_flip=4 if MIA327==2
replace MIA327_flip=2 if MIA327==3
replace MIA327_flip=1 if MIA327==4
replace MIA327_flip=3 if MIA327==5
tab MIA327 MIA327_flip
		*supervisor violates federal law, positive reward for complaintant and issue will be resolved
generate MIA328_flip=.
replace MIA328_flip=5 if MIA328==1
replace MIA328_flip=4 if MIA328==2
replace MIA328_flip=2 if MIA328==3
replace MIA328_flip=1 if MIA328==4
replace MIA328_flip=3 if MIA328==5
tab MIA328 MIA328_flip

	*clean DVs a second time: flip to run def not - def yes; exclude "don't know/not sure"
		*sexual harassment, female target
generate MIA320_flip_2=.
replace MIA320_flip_2=4 if MIA320==1
replace MIA320_flip_2=3 if MIA320==2
replace MIA320_flip_2=2 if MIA320==3
replace MIA320_flip_2=1 if MIA320==4
tab MIA320 MIA320_flip_2
		*sexual harassment, female target bystander
generate MIA321_flip_2=.
replace MIA321_flip_2=4 if MIA321==1
replace MIA321_flip_2=3 if MIA321==2
replace MIA321_flip_2=2 if MIA321==3
replace MIA321_flip_2=1 if MIA321==4
tab MIA321 MIA321_flip_2
		*sexual harassment, male target
generate MIA405_flip_2=.
replace MIA405_flip_2=4 if MIA405==1
replace MIA405_flip_2=3 if MIA405==2
replace MIA405_flip_2=2 if MIA405==3
replace MIA405_flip_2=1 if MIA405==4
tab MIA405 MIA405_flip_2
		*chemicals in water
generate MIA322_flip_2=.
replace MIA322_flip_2=4 if MIA322==1
replace MIA322_flip_2=3 if MIA322==2
replace MIA322_flip_2=2 if MIA322==3
replace MIA322_flip_2=1 if MIA322==4
tab MIA322 MIA322_flip_2
		*embezzling, bystander, $500-5k
generate MIA323_flip_2=.
replace MIA323_flip_2=4 if MIA323==1
replace MIA323_flip_2=3 if MIA323==2
replace MIA323_flip_2=2 if MIA323==3
replace MIA323_flip_2=1 if MIA323==4
tab MIA323 MIA323_flip_2
		*embezzling, bystander, $1k-50k
generate MIA407_flip_2=.
replace MIA407_flip_2=4 if MIA407==1
replace MIA407_flip_2=3 if MIA407==2
replace MIA407_flip_2=2 if MIA407==3
replace MIA407_flip_2=1 if MIA407==4
tab MIA407 MIA407_flip_2
		*embezzling, bystander, $2k-100k
generate MIA408_flip_2=.
replace MIA408_flip_2=4 if MIA408==1
replace MIA408_flip_2=3 if MIA408==2
replace MIA408_flip_2=2 if MIA408==3
replace MIA408_flip_2=1 if MIA408==4
tab MIA408 MIA408_flip_2
		*leaking
generate MIA324_flip_2=.
replace MIA324_flip_2=4 if MIA324==1
replace MIA324_flip_2=3 if MIA324==2
replace MIA324_flip_2=2 if MIA324==3
replace MIA324_flip_2=1 if MIA324==4
tab MIA324 MIA324_flip_2
		*exam cheating, rewards complaintant
generate MIA325_flip_2=.
replace MIA325_flip_2=4 if MIA325==1
replace MIA325_flip_2=3 if MIA325==2
replace MIA325_flip_2=2 if MIA325==3
replace MIA325_flip_2=1 if MIA325==4
tab MIA325 MIA325_flip_2
		*exam cheating, risk of losing job if reported
generate MIA326_flip_2=.
replace MIA326_flip_2=4 if MIA326==1
replace MIA326_flip_2=3 if MIA326==2
replace MIA326_flip_2=2 if MIA326==3
replace MIA326_flip_2=1 if MIA326==4
tab MIA326 MIA326_flip_2
		*supervisor threatens public interest, job security protected by law
generate MIA327_flip_2=.
replace MIA327_flip_2=4 if MIA327==1
replace MIA327_flip_2=3 if MIA327==2
replace MIA327_flip_2=2 if MIA327==3
replace MIA327_flip_2=1 if MIA327==4
tab MIA327 MIA327_flip_2
		*supervisor violates federal law, positive reward for complaintant and issue will be resolved
generate MIA328_flip_2=.
replace MIA328_flip_2=4 if MIA328==1
replace MIA328_flip_2=3 if MIA328==2
replace MIA328_flip_2=2 if MIA328==3
replace MIA328_flip_2=1 if MIA328==4
tab MIA328 MIA328_flip_2

*clean DVs a third time: exclude "don't know/not sure"; make 0 = "definitely not" or "probably not expose" and 1 = "definitely" or "probably expose"
		*sexual harassment, female target
generate MIA320_flip_3=.
replace MIA320_flip_3=1 if MIA320==1
replace MIA320_flip_3=1 if MIA320==2
replace MIA320_flip_3=0 if MIA320==3
replace MIA320_flip_3=0 if MIA320==4
tab MIA320 MIA320_flip_3
		*sexual harassment, female target bystander
generate MIA321_flip_3=.
replace MIA321_flip_3=1 if MIA321==1
replace MIA321_flip_3=1 if MIA321==2
replace MIA321_flip_3=0 if MIA321==3
replace MIA321_flip_3=0 if MIA321==4
tab MIA321 MIA321_flip_3
		*sexual harassment, male target
generate MIA405_flip_3=.
replace MIA405_flip_3=1 if MIA405==1
replace MIA405_flip_3=1 if MIA405==2
replace MIA405_flip_3=0 if MIA405==3
replace MIA405_flip_3=0 if MIA405==4
tab MIA405 MIA405_flip_3
		*chemicals in water
generate MIA322_flip_3=.
replace MIA322_flip_3=1 if MIA322==1
replace MIA322_flip_3=1 if MIA322==2
replace MIA322_flip_3=0 if MIA322==3
replace MIA322_flip_3=0 if MIA322==4
tab MIA322 MIA322_flip_3
		*embezzling, bystander, $500-5k
generate MIA323_flip_3=.
replace MIA323_flip_3=1 if MIA323==1
replace MIA323_flip_3=1 if MIA323==2
replace MIA323_flip_3=0 if MIA323==3
replace MIA323_flip_3=0 if MIA323==4
tab MIA323 MIA323_flip_3
		*embezzling, bystander, $1k-50k
generate MIA407_flip_3=.
replace MIA407_flip_3=1 if MIA407==1
replace MIA407_flip_3=1 if MIA407==2
replace MIA407_flip_3=0 if MIA407==3
replace MIA407_flip_3=0 if MIA407==4
tab MIA407 MIA407_flip_3
		*embezzling, bystander, $2k-100k
generate MIA408_flip_3=.
replace MIA408_flip_3=1 if MIA408==1
replace MIA408_flip_3=1 if MIA408==2
replace MIA408_flip_3=0 if MIA408==3
replace MIA408_flip_3=0 if MIA408==4
tab MIA408 MIA408_flip_3
		*leaking
generate MIA324_flip_3=.
replace MIA324_flip_3=1 if MIA324==1
replace MIA324_flip_3=1 if MIA324==2
replace MIA324_flip_3=0 if MIA324==3
replace MIA324_flip_3=0 if MIA324==4
tab MIA324 MIA324_flip_3
		*exam cheating, rewards complaintant
generate MIA325_flip_3=.
replace MIA325_flip_3=1 if MIA325==1
replace MIA325_flip_3=1 if MIA325==2
replace MIA325_flip_3=0 if MIA325==3
replace MIA325_flip_3=0 if MIA325==4
tab MIA325 MIA325_flip_3
		*exam cheating, risk of losing job if reported
generate MIA326_flip_3=.
replace MIA326_flip_3=1 if MIA326==1
replace MIA326_flip_3=1 if MIA326==2
replace MIA326_flip_3=0 if MIA326==3
replace MIA326_flip_3=0 if MIA326==4
tab MIA326 MIA326_flip_3
		*supervisor threatens public interest, job security protected by law
generate MIA327_flip_3=.
replace MIA327_flip_3=1 if MIA327==1
replace MIA327_flip_3=1 if MIA327==2
replace MIA327_flip_3=0 if MIA327==3
replace MIA327_flip_3=0 if MIA327==4
tab MIA327 MIA327_flip_3
		*supervisor violates federal law, positive reward for complaintant and issue will be resolved
generate MIA328_flip_3=.
replace MIA328_flip_3=1 if MIA328==1
replace MIA328_flip_3=1 if MIA328==2
replace MIA328_flip_3=0 if MIA328==3
replace MIA328_flip_3=0 if MIA328==4
tab MIA328 MIA328_flip_3

	*treatment effects
		*sexual harassment, female target
oneway MIA320_flip MIA320_treat, bonferroni
anova MIA320_flip MIA320_treat
margins MIA320_treat
			*post-hoc test for subject sex effects
anova MIA320_flip MIA320_treat##gender
			*FOR R3/response memo: replicate with OLS to show mathematical equivilance to ANOVA
reg MIA320_flip i.MIA320_treat
			*replicate w/ "don't know/not sure" excldued from analysis
anova MIA320_flip_2 MIA320_treat
			*replicate w/ 0/1 DV
anova MIA320_flip_3 MIA320_treat
		
		*sexual harassmwent, female target, bystander
oneway MIA321_flip MIA321_treat, bonferroni
anova MIA321_flip MIA321_treat
margins MIA321_treat
			*post-hoc test for subject sex effects
anova MIA321_flip MIA321_treat##gender
			*replicate w/ "don't know/not sure" excldued from analysis
anova MIA321_flip_2 MIA321_treat
			*replicate w/ 0/1 DV
anova MIA321_flip_3 MIA321_treat
		
		*harassment, male target (note: post-election)
oneway MIA405_flip MIA405_treat, bonferroni
anova MIA405_flip MIA405_treat
margins MIA405_treat
			*post-hoc test for subject sex effects
anova MIA405_flip MIA405_treat##gender
margins MIA405_treat gender MIA405_treat##gender
			*replicate w/ "don't know/not sure" excldued from analysis
anova MIA405_flip_2 MIA405_treat
			*replicate w/ 0/1 DV
anova MIA405_flip_3 MIA405_treat
		
		*chemicals in water
oneway MIA322_flip MIA322_treat, bonferroni
anova MIA322_flip MIA322_treat
margins MIA322_treat
			*replicate w/ "don't know/not sure" excldued from analysis
anova MIA322_flip_2 MIA322_treat
			*replicate w/ 0/1 DV
anova MIA322_flip_3 MIA322_treat
		
		*embezzling, bystander, $500-5k
oneway MIA323_flip MIA323_treat, bonferroni
anova MIA323_flip MIA323_treat
margins MIA323_treat
			*replicate w/ "don't know/not sure" excldued from analysis
anova MIA323_flip_2 MIA323_treat
			*replicate w/ 0/1 DV
anova MIA323_flip_3 MIA323_treat		
		
		*embezzling, bystander, $1k-50k
oneway MIA407_flip MIA407_treat, bonferroni
anova MIA407_flip MIA407_treat
margins MIA407_treat
			*replicate w/ "don't know/not sure" excldued from analysis
anova MIA407_flip_2 MIA407_treat
			*replicate w/ 0/1 DV
anova MIA407_flip_3 MIA407_treat
		
		*embezzling, bystander, $2k-100k
oneway MIA408_flip MIA408_treat, bonferroni
anova MIA408_flip MIA408_treat
margins MIA408_treat
			*replicate w/ "don't know/not sure" excldued from analysis
anova MIA408_flip_2 MIA408_treat
			*replicate w/ 0/1 DV
anova MIA408_flip_3 MIA408_treat
		
		*leaking re: immigration
oneway MIA324_flip MIA324_treat, bonferroni
anova MIA324_flip MIA324_treat
margins MIA324_treat
			*post-hoc test with pid3
				*make pid3 dem, rep, all others
generate dem_rep_oth=.
replace dem_rep_oth=1 if pid3==1
replace dem_rep_oth=2 if pid3==2
replace dem_rep_oth=3 if pid3==3
replace dem_rep_oth=3 if pid3==4
replace dem_rep_oth=3 if pid3==5
tab dem_rep_oth pid3
				*analysis
anova MIA324_flip MIA324_treat##dem_rep_oth
margins dem_rep_oth
pwcompare dem_rep_oth
			*replicate w/ "don't know/not sure" excldued from analysis
anova MIA324_flip_2 MIA324_treat
			*replicate w/ 0/1 DV
anova MIA324_flip_3 MIA324_treat
		
		*exam cheating 1
oneway MIA325_flip MIA325_treat, bonferroni
anova MIA325_flip MIA325_treat
margins MIA325_treat
			*replicate w/ "don't know/not sure" excldued from analysis
anova MIA325_flip_2 MIA325_treat
			*replicate w/ 0/1 DV
anova MIA325_flip_3 MIA325_treat
		
		*exam cheating 2
oneway MIA326_flip MIA326_treat, bonferroni
anova MIA326_flip MIA326_treat
margins MIA326_treat
			*replicate w/ "don't know/not sure" excldued from analysis
anova MIA326_flip_2 MIA326_treat
			*replicate w/ 0/1 DV
anova MIA326_flip_3 MIA326_treat
		
		*supervisor threatens public interest 1
oneway MIA327_flip MIA327_treat, bonferroni
anova MIA327_flip MIA327_treat
margins MIA327_treat
			*replicate w/ "don't know/not sure" excldued from analysis
anova MIA327_flip_2 MIA327_treat
			*replicate w/ 0/1 DV
anova MIA327_flip_3 MIA327_treat
		
		*supervisor threatens public interest 2
oneway MIA328_flip MIA328_treat, bonferroni
anova MIA328_flip MIA328_treat
margins MIA328_treat
			*replicate w/ "don't know/not sure" excldued from analysis
anova MIA328_flip_2 MIA328_treat
			*replicate w/ 0/1 DV
anova MIA328_flip_3 MIA328_treat

****
****


*within-subjects effects
	*clean DVs: flip to run, def not - def yes; put "don't know/not sure" in the middle
		*sexual harassment, female target
generate MIA320_flip=.
replace MIA320_flip=5 if MIA320==1
replace MIA320_flip=4 if MIA320==2
replace MIA320_flip=2 if MIA320==3
replace MIA320_flip=1 if MIA320==4
replace MIA320_flip=3 if MIA320==5
tab MIA320 MIA320_flip
		*sexual harassment, female target bystander
generate MIA321_flip=.
replace MIA321_flip=5 if MIA321==1
replace MIA321_flip=4 if MIA321==2
replace MIA321_flip=2 if MIA321==3
replace MIA321_flip=1 if MIA321==4
replace MIA321_flip=3 if MIA321==5
tab MIA321 MIA321_flip
		*sexual harassment, male target
generate MIA405_flip=.
replace MIA405_flip=5 if MIA405==1
replace MIA405_flip=4 if MIA405==2
replace MIA405_flip=2 if MIA405==3
replace MIA405_flip=1 if MIA405==4
replace MIA405_flip=3 if MIA405==5
tab MIA405 MIA405_flip
		*chemicals in water
generate MIA322_flip=.
replace MIA322_flip=5 if MIA322==1
replace MIA322_flip=4 if MIA322==2
replace MIA322_flip=2 if MIA322==3
replace MIA322_flip=1 if MIA322==4
replace MIA322_flip=3 if MIA322==5
tab MIA322 MIA322_flip
		*embezzling, bystander, $500-5k
generate MIA323_flip=.
replace MIA323_flip=5 if MIA323==1
replace MIA323_flip=4 if MIA323==2
replace MIA323_flip=2 if MIA323==3
replace MIA323_flip=1 if MIA323==4
replace MIA323_flip=3 if MIA323==5
tab MIA323 MIA323_flip
		*embezzling, bystander, $1k-50k
generate MIA407_flip=.
replace MIA407_flip=5 if MIA407==1
replace MIA407_flip=4 if MIA407==2
replace MIA407_flip=2 if MIA407==3
replace MIA407_flip=1 if MIA407==4
replace MIA407_flip=3 if MIA407==5
tab MIA407 MIA407_flip
		*embezzling, bystander, $2k-100k
generate MIA408_flip=.
replace MIA408_flip=5 if MIA408==1
replace MIA408_flip=4 if MIA408==2
replace MIA408_flip=2 if MIA408==3
replace MIA408_flip=1 if MIA408==4
replace MIA408_flip=3 if MIA408==5
tab MIA408 MIA408_flip
		*leaking
generate MIA324_flip=.
replace MIA324_flip=5 if MIA324==1
replace MIA324_flip=4 if MIA324==2
replace MIA324_flip=2 if MIA324==3
replace MIA324_flip=1 if MIA324==4
replace MIA324_flip=3 if MIA324==5
tab MIA324 MIA324_flip
		*exam cheating, rewards complaintant
generate MIA325_flip=.
replace MIA325_flip=5 if MIA325==1
replace MIA325_flip=4 if MIA325==2
replace MIA325_flip=2 if MIA325==3
replace MIA325_flip=1 if MIA325==4
replace MIA325_flip=3 if MIA325==5
tab MIA325 MIA325_flip
		*exam cheating, risk of losing job if reported
generate MIA326_flip=.
replace MIA326_flip=5 if MIA326==1
replace MIA326_flip=4 if MIA326==2
replace MIA326_flip=2 if MIA326==3
replace MIA326_flip=1 if MIA326==4
replace MIA326_flip=3 if MIA326==5
tab MIA326 MIA326_flip
		*supervisor threatens public interest, job security protected by law
generate MIA327_flip=.
replace MIA327_flip=5 if MIA327==1
replace MIA327_flip=4 if MIA327==2
replace MIA327_flip=2 if MIA327==3
replace MIA327_flip=1 if MIA327==4
replace MIA327_flip=3 if MIA327==5
tab MIA327 MIA327_flip
		*supervisor violates federal law, positive reward for complaintant and issue will be resolved
generate MIA328_flip=.
replace MIA328_flip=5 if MIA328==1
replace MIA328_flip=4 if MIA328==2
replace MIA328_flip=2 if MIA328==3
replace MIA328_flip=1 if MIA328==4
replace MIA328_flip=3 if MIA328==5
tab MIA328 MIA328_flip

	*figure 15: pooled means for 12 scenarios, regardless of experimental condition
		*sexual harassment, female target
means MIA320_flip
		*sexual harassmwent, female target, bystander
means MIA321_flip
		*harassment, male target (note: post-election)
means MIA405_flip
		*chemicals in water
means MIA322_flip
		*embezzling, $500-5k
means MIA323_flip
		*embezzling, $1k-50k (note: post-election)
means MIA407_flip
		*embezzling, $2k-100k (note: post-election)
means MIA408_flip
		*leaking
means MIA324_flip
		*fraud, retaliation, reform, reward
means MIA325_flip
		*fraud, retaliation, public interest
means MIA326_flip
		*severity of crime, job security
means MIA327_flip
		*crime, reward, punishment
means MIA328_flip

	*clean up names of the DVs for conversion to long form (here the rep meas are all 12 scenarios separately)
	*note: this renaming with the wb prefix and number after allows for the conversion to longform (need to have same prefix, and be numbered sequentially)
		*1 sexual harassment, female target
generate wb1=MIA320_flip
		*2 sexual harassmwent, female target, bystander
generate wb2=MIA321_flip
		*3 harassment, male target (note: post-election)
generate wb3=MIA405_flip
		*4 chemicals in water
generate wb4=MIA322_flip
		*5 embezzling, $500-5k
generate wb5=MIA323_flip
		*6 embezzling, $1k-50k (note: post-election)
generate wb6=MIA407_flip
		*7 embezzling, $2k-100k (note: post-election)
generate wb7=MIA408_flip
		*8 leaking
generate wb8=MIA324_flip
		*9 fraud, retaliation, reform, reward
generate wb9=MIA325_flip
		*10 fraud, retaliation, public interest
generate wb10=MIA326_flip
		*11 severity of crime, job security
generate wb11=MIA327_flip
		*12 crime, reward, punishment
generate wb12=MIA328_flip

	*generate long form data set for repeated measures anova
reshape long wb, i(caseid) j(scenario)
		*label new var scenario
label define scenario_label 1 "sh female" 2 "sh female, byst" 3 "sh male" 4 "h2o chems" 5 "embezz 1" 6 "embezz 2" 7 "embezz 3" 8 "leak" 9 "fraud 1" 10 "fraud 2" 11 "crime 1" 12 "crime 2"
label values scenario scenario_label
tab scenario

	*description of figure 15 tests in text of paper: repeated measures anova to test difference between the pooled means
set matsize 10000
anova wb caseid scenario, repeated(scenario)
margins scenario
pwcompare scenario

	*now post-hoc look at why reporting sexual harassment is lower than all other scenarios
		*get culled data set
clear
use "/Users/caseyklofstad/Dropbox/2018 CCES UM/paper 1 w-blowing and metoo/WB replication data.dta"
		
		*clean DVs: flip to run, def not - def yes; put "don't know/not sure" in the middle
			*sexual harassment, female target
generate MIA320_flip=.
replace MIA320_flip=5 if MIA320==1
replace MIA320_flip=4 if MIA320==2
replace MIA320_flip=2 if MIA320==3
replace MIA320_flip=1 if MIA320==4
replace MIA320_flip=3 if MIA320==5
tab MIA320 MIA320_flip
			*sexual harassment, female target bystander
generate MIA321_flip=.
replace MIA321_flip=5 if MIA321==1
replace MIA321_flip=4 if MIA321==2
replace MIA321_flip=2 if MIA321==3
replace MIA321_flip=1 if MIA321==4
replace MIA321_flip=3 if MIA321==5
tab MIA321 MIA321_flip
			*sexual harassment, male target
generate MIA405_flip=.
replace MIA405_flip=5 if MIA405==1
replace MIA405_flip=4 if MIA405==2
replace MIA405_flip=2 if MIA405==3
replace MIA405_flip=1 if MIA405==4
replace MIA405_flip=3 if MIA405==5
tab MIA405 MIA405_flip
			*chemicals in water
generate MIA322_flip=.
replace MIA322_flip=5 if MIA322==1
replace MIA322_flip=4 if MIA322==2
replace MIA322_flip=2 if MIA322==3
replace MIA322_flip=1 if MIA322==4
replace MIA322_flip=3 if MIA322==5
tab MIA322 MIA322_flip
			*embezzling, bystander, $500-5k
generate MIA323_flip=.
replace MIA323_flip=5 if MIA323==1
replace MIA323_flip=4 if MIA323==2
replace MIA323_flip=2 if MIA323==3
replace MIA323_flip=1 if MIA323==4
replace MIA323_flip=3 if MIA323==5
tab MIA323 MIA323_flip
			*embezzling, bystander, $1k-50k
generate MIA407_flip=.
replace MIA407_flip=5 if MIA407==1
replace MIA407_flip=4 if MIA407==2
replace MIA407_flip=2 if MIA407==3
replace MIA407_flip=1 if MIA407==4
replace MIA407_flip=3 if MIA407==5
tab MIA407 MIA407_flip
			*embezzling, bystander, $2k-100k
generate MIA408_flip=.
replace MIA408_flip=5 if MIA408==1
replace MIA408_flip=4 if MIA408==2
replace MIA408_flip=2 if MIA408==3
replace MIA408_flip=1 if MIA408==4
replace MIA408_flip=3 if MIA408==5
tab MIA408 MIA408_flip
			*leaking
generate MIA324_flip=.
replace MIA324_flip=5 if MIA324==1
replace MIA324_flip=4 if MIA324==2
replace MIA324_flip=2 if MIA324==3
replace MIA324_flip=1 if MIA324==4
replace MIA324_flip=3 if MIA324==5
tab MIA324 MIA324_flip
			*exam cheating, rewards complaintant
generate MIA325_flip=.
replace MIA325_flip=5 if MIA325==1
replace MIA325_flip=4 if MIA325==2
replace MIA325_flip=2 if MIA325==3
replace MIA325_flip=1 if MIA325==4
replace MIA325_flip=3 if MIA325==5
tab MIA325 MIA325_flip
			*exam cheating, risk of losing job if reported
generate MIA326_flip=.
replace MIA326_flip=5 if MIA326==1
replace MIA326_flip=4 if MIA326==2
replace MIA326_flip=2 if MIA326==3
replace MIA326_flip=1 if MIA326==4
replace MIA326_flip=3 if MIA326==5
tab MIA326 MIA326_flip
			*supervisor threatens public interest, job security protected by law
generate MIA327_flip=.
replace MIA327_flip=5 if MIA327==1
replace MIA327_flip=4 if MIA327==2
replace MIA327_flip=2 if MIA327==3
replace MIA327_flip=1 if MIA327==4
replace MIA327_flip=3 if MIA327==5
tab MIA327 MIA327_flip
			*supervisor violates federal law, positive reward for complaintant and issue will be resolved
generate MIA328_flip=.
replace MIA328_flip=5 if MIA328==1
replace MIA328_flip=4 if MIA328==2
replace MIA328_flip=2 if MIA328==3
replace MIA328_flip=1 if MIA328==4
replace MIA328_flip=3 if MIA328==5
tab MIA328 MIA328_flip
		
		*generate ave responses SH vignettes ("awb" = average whistle blow)
generate awb1=(MIA320_flip+MIA321_flip+MIA405_flip)/3
tab awb1
tabstat awb1, stats(mean semean)

		*generate ave repsonses for all other vignettes  ("awb" = average whistle blow)
generate awb2=(MIA322_flip+MIA323_flip+MIA407_flip+MIA408_flip+MIA324_flip+MIA325_flip+MIA326_flip+MIA327_flip+MIA328_flip)/9
tab awb2
tabstat awb2, stats(mean semean)

		*paired t-test for diff
ttest awb1==awb2

		*generate long form data set for repeated measures anova
reshape long awb, i(caseid) j(scenario)
			*label new var scenario
label define scenario_label 1 "sex harr" 2 "all others"
label values scenario scenario_label
tab scenario

		*rep meas anova with mixed for include gender (between), scenario (rep'd), interaction of the two, and caseid as control
anova awb gender / caseid|gender scenario scenario#gender, repeated(scenario)
margins scenario#gender, asbalanced emptycells(reweight)
pwcompare scenario#gender, asbalanced emptycells(reweight)

		*repeat above just DEMs
keep if pid3==1
anova awb gender / caseid|gender scenario scenario#gender, repeated(scenario)

		*rep above just REPs
keep if pid3==2
anova awb gender / caseid|gender scenario scenario#gender, repeated(scenario)







